



On the infrared behavior of Green's functions in Yang-Mills theory 
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Non-perturbative properties of QCD, such as color confinement, are encoded in the infrared behavior of 
correlation functions, e.g. propagators and vertices. Various analytic predictions have been suggested for these 
quantities in various gauges. Here we numerically test these predictions using lattice gauge theory. In particular, 
we present results for the 2- and 3-point functions for SU (2) Landau-gauge Yang-Mills theory in three and in 
four dimensions. Special attention is paid to systematic finite-volume effects. The gluon and ghost propagators 
are also evaluated in the so-called interpolating gauge (between the Landau and the Coulomb gauge), in order 
to study their gauge-dependence. Finally, we consider these propagators in Landau gauge at finite temperature, 
with the aim of understanding the effect of the deconfinement phase transition on their infrared behavior. All 
our results are compatible with the so-called Gribov-Zwanziger confinement scenario. 
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I. INTRODUCTION 

The confinement phenomenon in QCD is a long-range ef- 
fect and one can expect its explanation to be encoded in the 
infrared (IR) behavior of the QCD correlation functions. In- 
deed, various confinement scenarios make predictions for this 
behavior Q|. Unfortunately, a quantitative determination of 
these correlation functions is a very difficult task. 

Continuum methods, based e.g. on Dyson-Schwinger equa- 
tions (DSEs) Ol|3,S0L on the renormalization group j5[| or 
on effective potentials |6], lead to results that seem to agree 
with the Gribov-Zwanziger confinement scenario jj, 0] and 
with the Kugo-Ojima scenario |8]. These methods have the 
disadvantage of introducing approximations, which in general 
cannot be controlled and verified. On the other hand, numeri- 
cal studies using lattice gauge theory allow a determination of 
QCD correlation functions in the IR limit using approxima- 
tions that can in principle be controlled and quantified I9t ll0ll . 
In this case, since one has to work with finite volumes N , the 
system considered has a natural IR cutoff p ~ 1/ (aN), mak- 
ing difficult the study of the true IR limit. (This topic will 
be discussed in more detail in Section UTI) Clearly, a close 
cooperation between continuum and lattice methods can be 
a promising approach. In particular, lattice gauge theory al- 
lows to check the approximations usually considered in con- 
tinuum methods, or at least to provide constraints for possible 
Ansatze. The most important test is probably the verification 
of the approximations employed for the 3-point vertices when 
solving DSEs. Results for the three-gluon and the ghost-gluon 
vertices will be reported in SectionlHll 

Most of the studies of QCD correlations functions, both in 
the continuum and on the lattice, have been done in Landau 
gauge. On the other hand, since these quantities are gauge 
dependent, it is interesting to evaluate them in other gauges 
too (e.g. Coulomb gauge 1 11] and maximally Abelian gauge 
fl^l ). In principle, the explanation of the confinement mecha- 
nism based on the IR properties of the QCD correlation func- 
tion can also be gauge-dependent. Particularly important are 
the so-called interpolating gauges fl3l fl4ll . since they allow 



to relate different gauges and possible different confinement 
mechanisms. In Section fTVl we report 3d results for the gauge 
interpolating between the Landau and the Coulomb gauge. 

Finally, a clear understanding of the confinement mecha- 
nism could be useful in the study of the QCD phase diagram 
and of the deconfining phase transition 1 15]. In particular, one 
might expect to observe a modification of the IR properties 
of the correlation functions at finite temperature T. As a first 
step in this direction, the gluon and ghost propagators at finite 
T are considered (in the Ad Landau gauge) in SectionM 



H. FINITE-VOLUME EFFECTS 

Confinement is induced by long-range correlations. More 
precisely, the correlation length should be infinite 1 1]. Thus, 
when considering QCD correlation functions in momentum 
space, confinement should be related to their IR properties. 
This is the case for the Gribov-Zwanziger and the Kugo- 
Ojima confinement scenarios. Indeed, they predict a specific 
IR behavior for the gluon and the ghost propagators in Landau 
gauge: the gluon propagator should be IR suppressed — and 
it is expected to vanish at zero momentum — while the ghost 
propagator should be IR enhanced compared to the one of a 
massless particle. 

As said above, the IR-cutoff 1 / (aN) makes it difficult to 
study the IR properties of the correlation function and one 
should pay attention to possible finite-size effects. Of course, 
this analysis can be computationally very demanding. In this 
case, the consideration of three-dimensional lattices can be 
useful in order to access large lattice volumes. Let us note 
that the Gribov-Zwanziger scenario applies also to the three- 
dimensional case |4] and that results for the IR behavior of the 
propagators using DSEs are available for d = 3 ikl llfJl . 

In Figure [2 we report the comparison of DSE results with 
lattice data for the gluon and the ghost propagators. In the first 
case, finite-volume effects are evident for the smallest mo- 
menta. In particular, as the volume increases, the propagator 
becomes more and more IR suppressed. At the same time, the 
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FIG. 1 : Comparison of lattice data and DSE results for gluon D(p) (left) and ghost Dq (p) (right) propagators in the 3d case in Landau gauge 
as a function of the momentum p. The solid line is from DSE results 1 16]. Lattice simulations have been done at P = 4.2 considering various 
lattice volumes. Open circles correspond to lattice volume 20 3 [V ~ (3.5 fm) 3 ] U^l . open squares represent 30 3 [V ~ (5.2 fm) 3 ] flUl . open 
triangles are used for 40 3 [V ~ (6.9 fm) 3 ] fl^l . full circles correspond to 80 3 [V ~ (14 fm) 3 ] [20] and full triangles represent lattice volume 
140 3 [V « (24 fm) 3 ] HI. 



difference between two different lattices with a given volume 
ratio moves further and further into the IR limit with increas- 
ing volume. Thus, momenta of the order of 2sm(n/N)/a fss 
2%/{aN) are affected by finite-volume effects. From Figure^ 
one also sees that DSE and lattice results are in qualitatively 
agreement. Let us remark that DSE studies predict I17ll that 
the true IR regime is reached only below a certain momen- 
tum scale A/ < 200 MeV, i.e. on a finite volume one should 
have a sufficient number of non-zero momenta in the range 
2%/(aN) <C p < A/. With a lattice spacing a of order of 1 
GeV -1 one needs N 3> 50. Unfortunately, this regime is just 
reached by the lattices in the 3d case but not yet in Ad. 

The ghost propagator is much less affected by finite-volume 
effects, even though small effects can be found for the dress- 
ing function at small momenta I19I1 . Moreover, the compar- 
ison of the lattice data to DSE results is quite good (see Fig. 
[3- One should however recall that the IR enhancement of 
the ghost propagator seems to be related to an appropriate 
sampling of the so-called exceptional configurations Il0lll8tl . 
Thus, results obtained considering a relatively small statistics 
are likely underestimating this IR enhancement. 

We conclude that the smallest non-zero lattice momenta - 
of the order of 2%/ (aN) — are probably not reliable in lattice 
studies and should be treated with great care. A more thor- 
ough discussion of finite-size effects will be presented in 1 19]. 



IIL THREE-POINT VERTICES 

Besides the propagators, also the vertices are of great im- 
portance in understanding the IR properties of QCD. In par- 



ticular, assumptions on their IR behavior are crucial inputs 
for the analytic methods 0,0, 0]- At the same time, predic- 
tions for their IR behavior have been obtained using contin- 
uum methods 1 221 1231 12411 . As stressed in the Introduction, lat- 
tice gauge theory provides us with the possibility of verifying 
these results as well as the assumptions considered in analytic 
studies. The simplest of these vertices are the 3-point ones. 
In Landau-gauge pure Yang-Mills theory one can consider the 
ghost-gluon vertex and the three-gluon vertex. Kinematically 
these vertices are much more complicated than the propaga- 
tors, since they depend on three independent kinematic vari- 
ables. In addition, they have a complex tensor structure (espe- 
cially the three-gluon vertex). 

In order to simplify the presentation, we consider here only 
one specific kinematic configuration, namely the one in which 
two of the incoming momenta are orthogonal with respect to 
each other (see 1 1 8 , 25 ] for details and for other possible kine- 
matic configurations). Furthermore, instead of the various ten- 
sor elements we evaluate the quantity 11811 

r^l,abc/-'abe 
YH,abcj^adj^bejycfTlA,def' ^ ' 

Here the indices a,...,f are generic multi-indices, encom- 
passing field-type, Lorentz and color indices. Also, D ab are 
the propagators of the fields, G" bc represent the full Green's 
functions and r A abc are the corresponding tree-level vertices. 
This quantity is defined such that it becomes equal to one if 
the full and the tree-level vertex coincide. Furthermore, the ra- 
tio corresponds to the kernel of a one-loop not-amputated 
self-energy diagram in DSEs over the same kernel with full 
vertices replaced by the bare ones. Thus, the quantity G is a 
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FIG. 2: Top-left panel: ghost-gluon vertex at zero gluon momentum; open triangles are used for V = 30 3 1 25] and open circles for V = 30 3 (in 
both cases (3 = 4.2), full triangles are used for V = 20 3 1 18] and full circles for V = 30 3 (in both cases (3 = 6.0). Top-right panel: ghost-gluon 
vertex with the gluon and the ghost momenta orthogonal to each other; the graph interpolates data obtained from a 40 3 lattice at (3 = 4.2 1 25] . 
Again, the spike at zero momentum is an artifact, as the vertex function cannot be evaluated there. Bottom-left panel: 3-gluon vertex with 
one external momentum vanishing 1 18]; full symbols correspond to (3 = 4.2 and open symbols to (3 = 6.0; circles are used for V = 20 3 and 
triangles for V = 30 3 . Bottom-right panel: 3-gluon vertex with two external momenta orthogonal to each other 1 18]; the graph interpolates 
data obtained from a 30 3 lattice at (3 = 4.2 and at (3 = 6.0. Note that the spike at zero momentum is an artifact of the interpolating algorithm, 
since the vertex cannot be evaluated there. 



direct measure of how good the approximations in DSE cal- 
culations are when the bare vertices are used. 

Results for the quantity Q in three and in four dimensions 
for both vertices are shown in Figures |2] and [3] respectively. 
The ghost-gluon vertex is essentially constant for all momenta 



studied here as well as for other kinematic configurations 
I25II26 ! both in 3d and in Ad. In particular, there is no sign of 
enhancement or suppression even for the smallest momenta 
available. This is in a gre ement with results obtained using 
continuum methods I 23H24B . In addition, the constancy of this 
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FIG. 3: Same as Figure[2]but in the four-dimensional case. Here the lattice volume is 16 4 and simulations have been done at (3 = 2.2 12^1 . 



vertex confirms the most important assumption made in DSE 
studies, i.e. the use of a bare ghost-gluon vertex | Ill2l ll6ll24ll . 

The evaluation of three-gluon vertex is numerically much 
more demanding than that of the ghost-gluon vertex. In par- 
ticular, very large statistics are required for its investigation 
OH HH. especially in the three-dimensional case. Our results 
show that this vertex is suppressed at intermediate momenta, 
especially for d = 3. On the other hand, with our data, it is dif- 
ficult to understand completely its precise IR behavior. Let us 
note that, if this suppression were to be confirmed by studies 
using larger lattices, this result would contradict predictions 



obtained from DSE studies Nevertheless, it would 

not invalidate the assumptions usually made in the continuum 
calculations of the propagators Iulz [l6l.|24ll . 



IV. INTERPOLATING GAUGES 

All results presented so far have been obtained in the 
(minimal) Landau gauge. As said in the Introduction, it is 
eventually desirable to understand confinement in a gauge- 
independent way, if this is possible. Thus, the IR properties 
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Spatial gluon propagator D" | 




Temporal gluon propagator D 




FIG. 4: The spatial (top) and temporal (bottom) gluon propagator in 
the interpolating gauge J2j at X = 1/2 as a function of the temporal po 
and spatial \p\ momenta, from a 60 3 lattice at (3 = 4.2 [29]. Results 
have been interpolated; the spike at zero is an artifact of the missing 
point there. 



of the QCD correlation functions in other gauges are also in- 
teresting. Among other possible gauge-fixing conditions, the 
so-called interpolating gauges allow to interpolate smoothly 
between specific gauge conditions. One example is given by 
the gauge condition I13L l27ll 

XdoA b + dA b = 0. (2) 

This gauge interpolates between the Landau gauge (X = 1) and 
the Coulomb gauge (X = 0), even though the Coulomb gauge 
limit itself is likely not to be a smooth one [ 28 ] . Nonethe- 
less, this gauge allows one to study the Gribov-Zwanziger 
scenario when moving away from the Landau gauge. Since 
for X^l the <9(4)-invariance is explicitly broken, it is natural 
to consider two scalar gluon propagators l28ll . i.e. the time- 
time component D c q 'q/(N 2 — 1) and the 3-dimensional trans- 
verse one D Xr = D""/((d - l)(N 2 - 1)). In addition, the prop- 



| Ghost dressing function | 




FIG. 5: The ghost dressing function D(j(p)(pQ + p 2 ) in the inter- 
polating gauge J2j at X = 1/2 as a function of the temporal po an d 
spatial |p| momenta, from a 20 3 lattice at (3 = 4.2 |29T| . Numerical 
data have been interpolated. 



agators depend independently on the temporal momentum po 
and on the spatial momentum \p\. 

Results at X = 1/2 are shown in Figure for these two 
gluon propagators and in Figure |5] for the scalar ghost dress- 
ing function Dq{p){p 2 ) + p 2 ). Note that Doo(/?o,0) vanishes 
exactly for any po due to the gauge condition, even though 
(due to technical reasons) this is not clear from Figure|4] Both 
gluon propagators are found to be IR suppressed when con- 
sidered as a function of p\ + p 2 . Also, both propagators ex- 
hibit a clear maximum. Furthermore, both propagators are 
not suppressed as long as either po or \p\ are not sufficiently 
small. At the same time, the ghost propagator is IR enhanced. 
The difference in behavior at large pure-temporal and at large 
pure-spatial momenta is a direct consequence of how X ap- 
pears in the tree-level propagators. Thus, the IR behavior of 
these propagators is (qualitatively) the same as that found in 
Landau gauge. This is in agreement with results from DSE 
studies 1 28]. On the other hand, since the limit X — ► 1 is 
smooth, one wouldn't expect a strong deviation from the Lan- 
dau behavior for a value X = 1/2. Results for the non-smooth 
limit X — > will be presented elsewhere l2?jl . 



V. FINITE TEMPERATURE 

Yang-Mills theory undergoes a deconfining phase transition 
at finite temperature 1 15] and one expects this transition to 
change the confining properties of the theory. At the same 
time, if confinement is manifest in the IR properties of the 
QCD propagators, these might also change across the phase 
transition. At finite temperature, the gluon propagator can no 
longer be described by a single tensor structure and one needs 
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FIG. 6: The spatially transverse (top) and spatially longitudinal (bot- 
tom) gluon propagators at various temperatures 11911 . Full circles are 
at T — on a 20 4 lattice, full triangles are at T ~ 149 MeV on a 
10 x 26 3 lattice, open circles are at T « 298 MeV on a 4 x 34 3 lat- 
tice, and open triangles are at T ~ 597 MeV on a 2 x 42 3 lattice. All 
results have been obtained using (3 = 2.3. 



two independent structures 1 3 



j^ab pT p)Tab 



Lab 



Here, Pj N is purely spatial and transverse while P. 



P T 



i.e 



flV 

pL 



complements PL, to a four-dimensional-transverse 



<9(4)-invariant tensor-structure P pV = 8^ v — P/jPv/p 2 - In ad- 
dition, P^ v is longitudinal in the spatial sub-space. Further- 
more, all functions depend on the spatial momentum \p\ and 
on the energy po separately. The latter is discrete with values 
po = 2nTn, where n is an integer. Infrared or soft modes are 
those with zero energy, i.e. n = 0. 

Results for the soft modes of the two gluon propagators 
and of the (scalar) ghost propagator are shown in Figures [6] 
andQfor temperatures above and below the phase transition. 
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FIG. 7: Same as figure [fj] but for the ghost propagator. Also the 
symbols are the same except for full circles, which refer to a 32 4 
lattice at P = 2.3 (T = 0). 



Clearly, as expected, all propagators coincide with the zero- 
temperature perturbative tail for momenta p large compared 
to the temperature T. At the same time, each of the three 
propagators shows a different temperature dependence at low 
momenta. In particular, the 3af-transverse propagator D T (p) 
decreases smoothly (in the IR) with increasing temperature. 
At the same time, for the highest temperature considered here 
(T w 597 MeV), one sees a clear IR suppression with a maxi- 
mum for p w 600 MeV. Let us stress that, while at T — one 
does not see a decreasing propagator D T (p) even for a lat- 
tice volume of 52 4 l3"Tll . for this temperature above the phase 
transition the IR suppression is already manifest for a spatial 
volume of 42 3 . This suggests that, at high temperature, ei- 
ther the infrared suppression is stronger or that finite-volume 
effects are smaller, compared to the zero-temperature case. 
In the case of the 3ii-longitudinal propagator D L (p) we ob- 
serve an increase as the temperature is turned on starting from 
T = 0. However, the propagator decreases (in the IR region) 
when going from the second-highest to the highest tempera- 
ture. A possibility would be a propagator increasing (respec- 
tively decreasing) when the temperature rises for temperatures 
below (respectively above) the phase transition. Let us note 
that D L (p) does not seem to be IR suppressed, even though 
in the IR limit it does not grow as strongly as the propaga- 
tor of a massless particle 1 19]. Finally, the ghost propagator 
does not show any visible temperature dependence. Let us 
also stress that these results should be taken with care, es- 
pecially when considering the 3t/-longitudinal gluon propa- 
gator D L (p), since we find strong volume and discretization 
effectsCj. 

These results are in agreement with previous lattice l32ll 
and continuum results 13311 and they can be interpreted con- 
sidering a Gribov-Zwanziger-type scenario |19J. In particu- 
lar, they suggest that this confinement scenario prevails at all 
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temperatures, i.e. not only for small temperature, but also in 
the high-temperature phase. 

VI. SUMMARY 

Summarizing, we have presented results for gluon and 
ghost propagators considering the 3- and the 4-dimensional 
cases for various gauges at zero temperature. In all cases our 
results are in qualitatively agreement with the Kugo-Ojima 
and with the Gribov-Zwanziger scenarios. Moreover, an in- 
vestigation of the propagators in the interpolating gauge 
suggests that the Gribov-Zwanziger scenario is stable when 
moving away from the Landau gauge condition. 

At zero temperature, in Landau gauge, we have also eval- 
uated the ghost-gluon and the 3 -gluon vertices. This study 
provides additional hints and constraints for the assumptions 
usually considered in continuum studies. 



Finally, we have presented results for the two propagators 
at finite temperature around the deconfining phase transition. 
These results fit into an extension of the Gribov-Zwanziger 
scenario at finite temperature 

In all cases considered, the analysis of finite-volume effects 
should be done very carefully, especially for the gluon prop- 
agator. Thus, a complete understanding of the IR behavior 
of the QCD correlation functions and a confirmation (or dis- 
proof) of the available confinement scenarios will be obtained 
only when extending these studies to larger lattices. 
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